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ABSTRACT 

The  expected  wa i t  in  the  GI/G/1  queue  is  related  to  the  mean  and 
variance  of  the  idle  time.  For  arrival  distributions  which  are  IFR  or 
have  mean  residual  life  bounded  by  —,  simple  bounds  are  obtained  which 
give,  for  example,  the  expected  number  in  queue  to  within  at  most  one 
cus  tomer . 

By  equating  input  with  output,  relations  between  random  variables 
are  used  to  obtain  expressions  foi  the  moments  of  the  waiting  time  in 
terms  of  moments  of  the  inter-arrival,  service,  and  idle  time  distri¬ 
butions.  By  bounding  the  idle  time  moments,  bounds  are  obtained  on  the 
mean  and  variance  of  the  waiting  time,  the  mean  length  of  a  busy  period, 
and  the  probability  an  arrival  finds  the  system  empty.  Bounds  on  the 
mean  wait  lead  to  bounds  on  the  expected  virtual  wait. 

Similar  results  are  obtained  for  some  generalizations  of  the  GI/G/1 
queue,  including  batch  arrivals,  batch  service  and  priority  queues. 
Queues  where  the  first  customer  in  each  busy  period  has  some  added  delay 
are  also  considered. 

Some  preliminary  results  for  tandem  queues  are  given. 
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1  . 


Chapter  1 

SOME  RESULTS  AND  BOUNDS  FOR  ALL  GI/G/1  QUEUES 

1 .  I n  troduct i on 

Little  work  has  been  done  on  approximations  in  queue ing„  Emphasis  has 
been  on  complex  analytic  results.  Notable  exceptions  are  papers  by  Kingman 
(1962  (a)  and  (b))  and  recently  by  Newell  (1965).  The  work  of  Newell  is 
applied  primarily  to  traffic  light  problems,  whereas  Kingman's  is  more 
closely  related  to  this  thesis. 

In  this  chapter  some  new  results  are  found  for  various  indicators  of 
performance  in  the  GI/G/1  queue.  Bounds  which  are  easily  calculable  are  found 
for  such  items  as  the  expected  wa i t  in  queue,  expected  length  of  an  idle 
period  and  the  variance  of  inter-output  times. 

V/e  find  a  relationship  between  the  idle  time  between  busy  periods  and 
the  waiting  time  of  a  customer  in  queue.  The  expected  wait  in  queue  is  found 
in  terms  of  the  first  two  moments  of  the  inter-arrival,  service,  and  idle 
times.  For  Poisson  arrivals  the  idle  time  distribution  is  exponential,  and 
the  expected  wait  is  calculated  easily.  In  general  the  moments  of  the  idle 
distribution  are  difficult  to  calculate.  However,  an  upper  bound  for  all 
GI/G/1  queues  is  easily  found  in  terms  of  the  mean  and  variance  of  the  arrival 
and  service  streams  only  (see  also  Kingman,  1962  (a)  or  (b)).  A  lower 
bound  is  found  which  requires  knowledge  of  the  arrival  and  service  distri¬ 
butions,  and  not  just  the  first  two  moments. 

Only  stationary  queues  are  considered.  No  transient  results  are  given. 

2 .  Nota t i on 

The  following  notation  is  used  throughout  the  paper.  The  sign  ~  is 
used  to  signify  "with  distribution  function". 


2. 


V/e  shall  deal  exclusively  with  stationary  queues  in  this  paper  ,  by 
which  we  shall  mean  that  the  queueing  p-ocess  either  started  at  time  zero 
with  stationary  conditions  or  that  it  started  with  some  initial  condition 
(such  as  the  wait  in  queue  of  the  first  customer  is  zero)  but  that  time 
was  at  . 

By  the  subscript  n  (e.g.,  W  )  we  shall  be  referring  to  the  n-th  cus¬ 
tomers  in  a  stationary  stream.  V/hen  it  is  not  required  to  note  the  order 
of  the  customers  the  subscript  will  be  dropped. 


T 

n 


S 

n 


U 

n 


T 

n 


P 


W 

n 


V 


I 


B 


=  time  between  n-th  and  (n+l)-th  arrival,  Tn  ~A(t),  E[T  ]  =  1A 

=  service  time  of  n-th  customer  >  S  r^j  G(t),  E[S  ]  =  1/p, 

-  S  -  T  ,  U  ~  K(t) 
n  n  n  '  ' 

=  time  between  n-th  and  (n-i-l)-th  departure 
=  X/p , 

=  wait  in  queue  of  n-th  customer,  W  —  V/ ( t ) 

=  virtual  wait  in  a  stationary  queue  (see  section  4) 

=  length  of  idle  period  between  busy  periods,  I  ~  H(t) 

=  length  of  busy  period  ,  B  ~  B(t) 


It  is  possible  in  some  queueing  situations  that  an  arrival  and  service  can 
take  place  together,  leading  to  problems  in  defining  what  is  an  idle  period 
for  the  queue.  V/e  shall  define  P[l=0]  =  0,  and  thus  if  an  arrival  occurs 
at  the  instant  the  last  customer  present  departs,  the  busy  period  continues, 
and  ends  only  when  the  facility  is  empty  for  a  positive  length  of  time. 


=  number  served  in  a  busy  period 


3» 


D  =  total  delay  in  system  -  V/  +  S 
N  =  number  in  the  queue  at  a  random  point  in  time 

q 

=  n-th  moment  about  origin  of  random  variable  with  distribution  F. 

The  superscript  is  dropped  for  n=l  ,  e.g.,  =  1A,  =  1/p,  =  E[|] 

2 

=  variance  of  a  random  variable  with  distribution  F 
2  2  2 

Cf  =  or^/(v^)  ,  where  c^.  is  the  coefficient  of  variation 
=  P[Ar rival  finds  the  system  empty] 
pc(t)  =  1  -  F (t)  for  any  distribution  F 


3 •  The  Wait  in  Queue  and  the  Idle  Period 

By  equating  input  times  with  output  times  relationships  between  the 
moments  of  the  arrival,  service,  idle,  and  waiting  time  distributions  are 
now  found. 

The  following  result  is  given  in  Riordan  ( 1 962 ) ,  but  because  of  its 
usefulness  it  is  proved  here. 

Theorem  1,1;  For  all  GI/G/1  queues  with  p  <  1  , 


a  A 

o  h 


0) 


Proof :  Consider  the  time  between  the  n-th  and  (n-M)-th  departures 


T 

n 


H- 


D 


n+1 


=  S 


n+1 


+  X 

n 
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4. 


where 


X  =  0  i  f  W  ,  ,  >  0 
n  n+ 1 


-  -(D  -T  )  if  W  4l  =  0. 
'  n  n  n+1 


If  U  =  S  -T  we  have  the  fundamental  relationships 
n  n  n 


W  J  .  =  W  +  X  +  U 
n+1  n  n  n 


=  Max  [0,  W  +  U  ]. 
L  n  nJ 


When  Xn  >  0 ,  X^=l .  Hence,  taking  expectations  in  (2)  and  assuming 


stat  ionar  i  ty 


E[Xnl  =  aQE[ 1 ]  =  -E[Up]  which  proves  (1). 


As  examples,  for  Poisson  arrivals  aQ  =  (1-p)  and  the  idle  distribution 
is  exponential  with  mean  For  the  constant  arrival,  constant  service 
case  (D/D/1  queue)  a^  =  1  and  1  =  ^ 


The  following  expression  is  now  derived  for  the  expected  wa  i  t  in 


queue , 


Theorem  1.2;  For  all  Gl/G/1  queues  with  p  <  1, 


Ervn  -  E[»2]  .  iLlil 

tLWJ  ■  -2E[U'J  2E[I] 


2  ^  2,  ,2  (2) 
X  +  aq)  +  0_P)  vh 

2X{  1-p)  "  2^~ 


(4) 


5. 


Proof:  Write  (2)  as  W  ,  ,  -  X  =  W  +  U  .  Square  both  sides  and  note  that 
-  v  '  n+1  n  n  n 

Wn-HXn  “  °’  3ivin9 


2  2  2  2 
W  ,  .  +  X  =  W  +  2W  U  +  U  . 
n-H  n  n  n  n  n 


Taking  expectations,  since  W  and  Un  are  independent,  and 


E[XZ]  =  a  E[|2]  we  have 
n  o  J 


aoE[|2]  =  2E[Un]E[Wn0  +  E[U2]. 


Using  theorem  1.1  the  result  follows. 


It  is  important  to  note  the  special  v/ay  in  which  the  moments  of  the  idle 
Ef  |21 

distribution  occur.  — iS  the  mean  °f  ar>  equilibrium  excess  idle  distri¬ 
bution;  that  is,  it  is  the  mean  of  a  random  variable  with  distribution 

function  — du.  This  is  a  we  1 1  known  result  in  renewal  theory. 

°  "h  (2) 

v  , 

Consider  again  our  two  examples.  For  Poisson  arrivals  ^ 

o  Vh 

P(l+Cg) 

In  this  case  (4)  reduces  to  E[W’J  -  p)  ’  3  we  ^  known  result.  For  the 

?  ?  (?)  7112 

D/D/1  queue  cr^  =  =  0  and  \v  =  (v^)  =  (^-  -  — )  ,  in  v/hich  case  (4) 

reduces  to  E[W]  =  0. 

An  expression  for  the  variance  of  the  wa i t  is  now  found  in  a  similar 
manner  and  is  given  by 


Theorem  1.3:  For  all  Gl/G/1  queues  with  p  <  1, 


(5) 


T 


rf 


,3-, 


2-, 


2  Ef  u  1  f  Ej  U__l 

CTw  '  -3E[  UJ  "  L  -  2E  [  U] 


3-, 


2,-,  ?, 


+  ,'lLLli  _  r§u__ 

-3E[1]  L 2E[  I  jJ 


or 


°w 


X(v^-v^)  -i-  3  (ov^-v^)  rX2(o2-!-cr2)  -i-  (l-p)2n2 

v  q  a  '  ■’  v  a  q  '  ^  | _ __a _ g/  v  ^ 


where 


3(l-p) 


2  (2)  2 

C  u  “  V  /-v  , 

e  ,h  e  ,h  e  ,h 


v 


(3) 


2X(J~p) 
■<2>  2 


a 


e  ,h 


3v, 


rits _ n 

L2v,  J  * 


Proof:  V/rite  (2)  as  V/  ,  .  -  X  =  V/  +  U  and  cube  both  sides.  Note  that 
■ - —  n+1  n  n  n 

XVH  =  V/n-!-  lXn  =  0.  Using  (1)  and  (4)  after  taking  expectations 
the  resul t  fol lows . 

The  expression  for  the  expected  wait  is  of  particular  interest  in  queue¬ 
ing  and  it  is  seen  to  depend  only  on  the  first  two  moments  of  the  inter- 
arrival,  service,  and  idle  distributions.  In  general  these  idle  period  moments 
are  difficult  to  calculate  but  bounds  will  be  obtained  for  them  and  this  is 
the  subject  of  section  (7)  of  this  chapter  and  sections  ( 1  -X| )  of  Chapter  2. 


4 .  The  Virtual  V/a  i  t 

The  virtual  wait  is  defined  here  to  be  the. time  an  arrival  would  have 
to  v/a  i  t  in  queue  if  he  arrived  at  some  random  point  in  time  into  a  station¬ 
ary  queue.  The  expected  value  of  V  is  found  easily  in  terms  of  the  expected 
value  of  V/,  the  actual  wait  in  queue. 

First  we  show  that  at  a  random  point  in  time  in  a  stationary  queue, 
the  probability  the  system  is  empty  is  (1-p).  The  times  between  the  starts 
of  busy  periods  are  independent  and  identically  distributed,  say  with 
distribution  B  (t),  the  distribution  of  (BH).  Hence,  the  queue  has  an 
imbedded  renewal  process,  and  note  that  following  the  end  of  an  idle  period, 


y*.  »> 


7. 


the  distribution  of  the  next  busy  period  is  independent,  and  distributed 
as  B(t).  Let  M'(t)  be  the  expected  number  of  busy  cycles  (i.e.,  from  the 
start  of  one  busy  period  to  the  start  of  the  next)  in  (0,t]  starting  at 
t=0  with  the  start  of  a  busy  period. 


r  r 

P[Busy  at  Tl  =  J  B  (T-u)dM'(u) 


and  applying  the  Key  Renewal  theorem,  if  B*  (t)  is  non-lattice, 


1  ini  P[Busy  at  T]  =  F[B^+E['i']  =  p 


In  the  lattice  case,  the  Cesaro  limit  may  be  found  by  Tauberian  arguments. 

Now  let  V  =  X-t-Y,  where  X  is  the  excess  service  time  of  the  customer 
in  service  and  Y  is  the  sum  of  the  service  times  of  all  those  in  queue  when 
a  random  arrival  occurs.  Then 


E[N  ] 

E[Y]  =— pE[W]. 

P 


Conditioning  only  on  whether  or  not  a  random  arrival  finds  the  system  busy, 

(2) 

vv  1 

E[X]  :■ 


2v  ‘ 

g 


Hence , 


E[V]  = 


p(c;  +  o 


2ii 


+  pE[V/] 


(6) 


It  is  interesting  to  note  from  (6)  that  E[V"J  =  E[V/'J  if  and  only  if 


E[W]  = 


P  (C^  +  1) 
2P  ( 1 -P ) 


which  is  the  case  for  Poisson  arrivals.  Using  equation 


8. 

i 

(^i)  in  (6)  gives 


E[V'J 


7  7  7 

p(C  +  C  )  M-C 
r  q  a  a 

2n, (1-p)  _  +  "V 


2"h  ' 


Using  either  (6)  or  (7)  in  the  D/D/1  queue 


gives  the  Cesaro  mean  E[V] 


o 

V* 


5  •  The  Variance  of  the  Output 

It  is  obvious  that  E[t  1  -  E[T  ]  =  The  variance  is  found  as 

n J  n  X 

follows.  From  the  relationships  in  section  3  > 


Var  [t  ]  =  Var  [S  .]  +  Var  [X  ] 
nJ  n+lJ  nJ 

(G) 

and 

Var  [V/  ,  ,-X  ]  =  Var  [D  -T  1  =  a2  +  a2  +  o2 . 
n+1  nJ  nnJ  a  g  w 

(8a) 

But 

Var  [W  ,  ,  -X  1  =  o2  +  Var  [X]  -  2Cov  (W  ,  ,X  )  . 

L  n+1  nJ  w  L  J  '  n+1  n 

(8b) 

Now  V/ 

.  ,X  =  0  and  hence , 
n+1  n 

Cov  (W  .X  )  =  -E[V/](~  -  -)  . 
v  n+1  n  J  X  p, 

Using 

this  with  (8),  (8a),  and  (8b)  gives 

Var  [t  ]  =  o2  +  2o2  -  7-  (l-p)E[\7]. 

L  nJ  a  g  X  p  J 

(9) 

Using  equation  (h)  for  E[V/]  we  have  finally, 


Var  (t  ) 
'  n' 


~  a 


ibP-t 

x2 


+ 


Ihfil 

X 


(2) 


v. 


JL'U.Jl*-* 


-  * 


?  1 

For  the  M/G/1  queue  this  gives  Var  (t  )  =  aj  +  — ~£~,  so  that  the  variance 

n  9  X 

of  the  output  of  the  M/G/1  queue  is  known  exactly  when  the  mean  and  variance 


of  the  service  distribution  are  given. 

2  1  1 

If  G(t)  is  exponential,  a  =  ~  anc^  Var  (t  )  - 

9  X 


In  the  case  of 


constant  arrivals,  constant  service,  a  ~  0  and  Var  (t  )  =  0. 

g  n' 


6 .  The  Covariance  of  Consecutive  Outputs 

V/e  now  derive  an  expression  for  the  covariance  between  two  consecutive 
output  times.  This  expression  depends  on  the  covariance  between  a  service 
time  and  the  wait  of  the  next  customer.  However  bounds  can  be  obtained  to 
show  that  it  is  bounded  close  to  zero.  This  is  done  in  section  7. 

Theorem  1  A:  For  all  GI/G/1  queues  with  p  <  1, 


Cov  (t  ,t  .)  -  Cov  (S  .  ,W  l0)+r-(l-o)  -  , 

N  n »  n  i  l '  '  n+1  n+2  X  g 


n’  n-i- 


(10) 


r  i  r 

where  ct  =  Ej^Max  [ 0 , U n  1  -  x  dK(x) 


P  roof : 


t  =  T  -  D  +  D  ,  , 
n  n  n  n-M 


(t  T  ,  i )  ~  (T  -D  -I-D  , )  (T  ,  ,-D  ,  ,+D  l0)  , 
n  n-i-r  N  n  n  n-i-l/N  n-M  n-M  n-i-2' 


and  in  a  straight-forward  but  tedious  manner  we  find  that 


Cov  (t  ,t  -  E[  (D  -T  )D  ,  ,1  -  ET  (D  -T  )D  -  E[D  ,  ,]  +  F.[D  ,  ,D  .  (a) 

'  n  n+r  n  rr  n-M J  n  n'  n-l-2J  nil  n+1  n-i-2J  v  ' 


S  i nee  D  .  =  W  ,  +  S  ,  ,  and  W  ,  =  Max  [ 0  ,D  -T  ]  , 
n+1  n+1  n-H  n+1  n  nJ 


E[cL,] 


“  E[wn+d  2Et'-'n+|Sn+,l  +  (^)2- 


Using  (b)  and  (c)  in  (a),  after  some  cancellation, 


Cov  (t  ,t  =  E[(W  ,+S  ,  -  (D  -T  ))W  ^  -  a2. 

'  n  n+1  'n+1  n+1  '  n  n"  n-i-2 J  u,  g 


=  -E[(Hin(0,Dn-Tn))Wn+2]  +  Cov(Sn+, ,Wnt2)  -  a 


W  Min  (0  ,D  -T  )  =  Max  (0  ,W  +S  ,,-T  ..)  Min  (0,D  -T  ). 
n+2  '  n  n'  n+1  n+1  n+1  n  n' 


When  D  -T  <0,  then  W  .  =  0, 
n  n  n+1 


•••  Wn+2  Min  <°'W  =  [°>Sn+rW  Mi"  <°VTn> 


which  by  independence  gives 


E[Wr+2  Min  (0,Dn-Tn)]  »  -«E[Max  (0,Dp-Tn)  -  (Dn-Tp)] 


,1  K 
-  $ 


Putting  this  in  (d)  gives 


Cov(Tn)Tn+l)  =  Cov(Sn+i  »Wn+?)  H'  r  (1"o) 


n+ 1  ’  n+2/  X 


2 

‘  V 


11. 


which  establishes  (10). 


7.  Some  Bounds  for  All  GI/G/1  Queues 

Using  the  results  of  the  previous  sections  some  simple  bounds  can  be 
found  for  various  factors  in  the  GI/G/1  queue,  such  as  the  mean  length  of 
an  idle  period  and  the  mean  wait  in  queue, 

a)  The  mean  idle  time.  Since  aQ  g  1  ,  theorem  1.1  gives  a  lower  bound  on 
the  length  of  an  idle  period, 


(id 


(12) 


This  upper  bound  for  all  Gl/G/1  queues  is  also  found  by  Kingman 
(1962  (a)  and  (b)).  Equality  holds  for  the  D/D/1  queue. 

The  importance  of  these  bounds  is  that  they  involve  at  most  only 
the  first  two  moments  of  the  arrival  and  service  distributions  and  further 
knowledge  of  the  distributions  is  not  required.  However,  if  K(t)  is 
known  (or  alternatively  if  A(t)  and  G(t)  are  known)  a  lower  bound  on  the 
wa  i  t  in  queue  can  be  found  as  follows. 


! 


12, 


Theorem  1.5:  Let  t  be  a  solution  of 


x  =  1  KC(u)du,  x  5  0,  where  (Sn~Tn)  ~  K(t) 

which  exists  and  is  unique  if  and  only  if  p  <  1.  Then  for  all 
GI/G/1  queues,  E[W]  §  l. 

Proof :  Recall  the  fundamental  equation  (3) 


W  L ,  =  Max  [0,  W  +  U  ]. 
n+1  n  n 


Then  [Wn+1|Wn=x]  =  [Hex  0,  x+Uj  , 


and 


CO 

E[W  +i |W  =X1  =J  KC(u)du  all  x  ^  0. 
!  -x 


(13) 


Now  let 


K  (u)du  =  g(x)  which  is  a  continuous  convex  function  for 


-x 


x  5  0,  with  g'(x)  =  Kc(-x)  ,  so  Kc(0+)  =  g '  (0“)  =  P[Up  >  0]  and 
g 1  (x)  ->  1  as  x  ->  « . 

f° 

Let-0  -  E[Min  (0,U  )]  -  K(u)du  and 

n  J 

-CO 

CO 

o'  =  E[Max  (0,U  )]  «  J  KC(u)du. 


Then 


From  (13) 


c  1  1 

a  -  6  - - 

(t  X 


E[W„n'J  =  I0.9(x)d  Wnw 


rrzr-}*" 


pe 
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or  E[Wn+]]  =  E[g(V/n)]. 

Using  Jensen's  inequality  for  the  expected  value  of  a  convex 
function  of  a  non-negative  random  variable, 


E[Wn+1]  §  g(E[Wn]),  so  that 


E[W]  5  KC(u)du, 

J -E[W] 


Consider  the  equation 


(14) 


This  can  be  wr i tten 


r° 

x  =  a  +  J  KC(u)du  x  =  0. 

-x 


x  =  J  K°(u)  du,  x  ^  0 
-x 


The  situation  is  drawn  in  Figure  1.  The  equation  has  a  solution  if 
and  only  if  the  two  curves  cross.  If  a  =  0,  x=0  is  a  solution; 
if  a  >  0  the  curves  cross  if  and  only  if  for  x  sufficiently  large, 


x  >  cr  -f- 


r° 

(u)  du 


<=> 


r° 

K(u)du  >  cr 


or  i  f  and  only  i  f  (3  >  o'. 


-x  -x 

But  p  >  cr  if  and  only  if  r->  — .  Uniqueness  comes  from  convexity 

A  (X 

arguments.  Uniqueness  fails  only  when  the  two  curves  coincide  over 
some  range,  [a,b]  say.  This  implies  g 1  (x)  -=  K  (-x)  =  1  on 
[a,b]  ->  g'(x)  -  1  on  [a,®)  -->  curves  don't  cross.  In  the  case 
p  -  1  ,  either  no  solution  exists  or,  for  example  in  the  case  of  the 


D/D/1  queue,  an  infinite  number  of  solutions  exist  with  p  =  l  . 

So  for  p  <  1,  let  t  be  the  unique  solution  of  (14) .  It  is 
now  shown  that  l  ~  E[V/*J.  This  is  obvious  from  Figure  1  and  equations 
(13)  and  (14).  If  1=0  the  inequality  is  trivial.  If  t  >  0  then 

r° 

o'  >  0  and  for  all  0  5  x  <  t  ,  x  <  a  -I-  K  (u)du  from  the  uniqueness 

-x 

Pro 

property  of  t.  Hence,  if  E[W]  <  l  ,  then  E[W]  <  KC(u)du  which 

J -E[W] 

contradicts  (14)  and  the  theorem  is  proved. 

Summarizing,  we  have  shown  that  for  all  GI/G/1  queues  wi th  p  <  1 


15. 


2  2 

where  t  is  the  unique  solution  of  (14).  For  °a  *1'  °g  >  0  (i.e.,  all 

except' the  D/D/1  queue),  both  bounds  tend  to  infinity  as  --  -*  —  >  0. 

A.  pi 

However,  their  ratio  may  diverge  in  a  particular  case  as  is  shown 

below  for  the  case  of  M/M/1  queue. 

» 

For  the  Poisson  arrival,  Exponential  service  queue  it  is  found 

that 

k(t)  =  ^re"^  1 8  0 

=  tso 

which  gives 

KC(t)  =  t  S  0 


1  2 

l  -  - (—)  log^  (l-p  )  which  -»  °°  as  p  -*  1”. 

However,  it  is  easy  to  show  that  lim  (l-p)  log  - y  =  0  and 

p~r  e 

hence,  the  bounds  diverge.  The  upper  and  lower  bounds  and  true  value 
of  E[VfJ  are  shown  in  Figure  2  for  fixed  A-l  and  varying  p. 

c)  The  variance  of  the  output. 

The  variance  of  the  output  distribution  is  given  in  equation  (9). 
Using  arguments  similar  to  those  in  b)  the  following  upper  and  lower 


17. 


Hence,  Cov  (Sn+,Wn+2)  5  E[S^,]  -  ^ 

and  using  this  in  (10)  gives 

Cov  (t  ,t  . )  §  r-  ( 1  -p )  (cv  -  — ) . 

'  n  n+ r  A  '  p, ' 

It  is  easy  to  show  that  of  <  ~  and  thus  that  this  expression  is  negative. 
An  upper  bound  is  found  as  follows: 

Max  [0,W  ,,+S  Ll-T  «W  .,  +  Max  [0,S  Ul-T 

L  J  n-i-1  nil  n+lJ  n-H  n-H  n+1 J 

and  using  this  gives 

Cov  3  E[Sn+l  tO,Sn+I-Tn+,]]  S  E[Sn2+)], 


which  leads  to: 

Theorem  1,6;  For  all  GI/G/1  queues  with  p  <  1 


h.  1 


V> 


£  Cov 


(VTn+l> 


Of(r 


^  +T 


(17) 


The  lower  bound  is  negative  and  upper  bound  positive  so  in  general  no 

conclusion  can  be  drawn  as  to  the  sign  of  the  covariance.  However,  it 

has  been  bounded  to  within  r— . 

Xp 

e )  The  Virtual  Wait . 

The  bounds  obtained  in  part  (b)  and  given  in  equation  (15)  can  be 
used  in  equation  (6)  to  show  that,  for  all  GI/G/1  queues  with  p  <  1, 


2^ 


l]  S  E[V] 


p(C2  +  c2) 

-  q  a 

=  2n(l-p) 


18. 


(IB) 


r-v^KTYv-v-  —  •:*.  «-»* war*"*"  -  •*Mi-gs55r*****.n.«.'« 


19. 


r;cw»cr  rmajr,- 


Chapter  2 

QUEUES  WITH  ARRIVAL  DISTRIBUTIONS  WHICH  HAVE  MONOTONE  FAILURE  RATES 


1 .  I  ntroduct i on 

In  the  previous  chapter  it  was  seen  that  the  moments  of  the  idle  time 
distribution  occurred  in  many  of  the  expressions.  The  idle  time  distribution 
is  some  complicated  tail  distribution  of  an  inter-arrival  time  and  it  might 
be  conjectured  that  by  placing  some  restriction  on  the  inter-arrival  time 
distribution  one  might  obtain  some  desirable  properties  of  the  oments  of 
the  idle  period  regardless  of  the  service  distribution.  This  indeed  turns 
out  to  be  true. 

The  first  restriction  to  be  placed  on  the  arrival  distribution  is  to 
restrict  it  to  the  class  with  decreasing  mean  residual  life  (DMRL) .  In  this 
chapter  the  words  decreasing  and  increasing  are  used  in  the  weak  sense  and 
always  should  be  read  to  mean  non- increas ing  and  non-decreasing  respectively. 
The  symbols  i  and  t  will  be  used  respectively  for  decreasing  and  increasing 
in  this  weak  sense. 

Def i n i t ion  1 .  A  non-discrete  distribution  F  has  DMRL  (IMRL)  if  and  only  if 


FC(x)dx  1 
FC(t)  (I) 


for  all  t  I  0  when  finite. 


The  expressions  and  symbols  in  parenthesis  should  be  read 
together . 

A  slightly  stronger  assumption  ori  the  arrival  distribution  will  also 
be  used  which  implies  the  above  assumption. 

Definition  2 .  A  non-discrete  distribution  F  has  increasing  failure  rate 
(is  I FR )  if  and  only  if 


F(M-A)  -  F(t)  1 
FC(0 
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for  A  >  0  and  all  t  ?  0  where  finite. 


If  F  is 


S  P 
n-k 


n 


d i screte  , 


then  it  is  1FR  if  and  only  if 


all  k  5  0  where  finite. 


It  is  said  to  be  DFR  (have  decreasing  failure  rate)  when  t 
is  replaced  by  l  in  the  above  expressions. 

These  concepts  are  widely  used  in  reliability  theory  where  strong 
physical  justifications  can  be  given  for  their  use  in  particular  problems. 

In  queueing  an  IFR  arrival  distribution  would  have  the  following  physical 
interpretation.  Given  it  has  been  a  time  t  since  the  last  customer  arrived, 
the  probability  that  a  customer  arrives  in  the  next  small  interval  A  is 
increasing  in  t.  Besides  any  physical  justification  many  parametric 
families  have  this  property;  for  example  the  gamma  and  V/eibul  distributions 
in  certain  parameter  ranges,  and  the  truncated  normal  and  modified  extreme 
value  distributions.  The  degenerate  distribution  of  the  constant  arrival 
queue  also  has  the  IFR  property.  For  a  fuller  discussion  on  these  properties 
the  reader  should  consult  Chapter  2  of  R .  E.  Barlow  and  F.  Proschan  (1965). 

For  IFR/G/1  queues,  (that  is,  the  class  of  GI/G/1  queues  whose  arrival 
distributions  have  the  IFR  property)  it  is  shown  that  simple  expressions 
can  be  obtained  to  bound,  for  example,  the  expected  number  in  the  queue  to 
within  at  most  one  customer.  These  bounds  involve  only  the  mean  and 
variance  of  the  arrival  and  service  streams.  For  the  special  class  of 
D/G/l  ,  (constant  arrival,  general  service),  the  expected  number  in  the 
queue  is  bound  to  within  at  most  one  half. 


21. 


2 .  Some  Properties  of  the  Idle  Distribution 

In  this  section  two  theorems  are  proved,  which  give  some  of  the  useful 
properties  of  the  idle  distribution.  These  are  then  used  in  section  3  to 
obtain  bounds  for  certain  factors  in  the  given  class  of  queues.  In  what 
follows,  symbols  and  expressions  in  parenthesis  should  be  read  together. 

Theorem  2.1;  For  all  Gl/G/1  queues  where  A(t)  is  restricted  to  have  D M2 L 

(imrl)  , 

r"  ulixife  *  r-  A^cidx  alltB0. 

Jt  HC(t)  Jt  AC(t) 


The  inequalities  are  tight  for  the  exponential  distribution. 


Proof:  Let  X  -  time  from  the  last  arrival  to  enter  a  busy  period  to  the 
end  of  the  busy  period  and  let  X  ~  § (x)  (see  Figure  3). 
Note  that  by  definition  (see  section  1,  Chapter  1)  ?  (0  )  =  0. 


c 

E 

1-  (U 
QJ  ■>-’ 
jo  in 
E  >• 
3  in 


Figure  3: 


f- - T - - - 

A  Realization  Showing  the  End  of  One  Busy  Period  and 

the  Start  of  the  Next. 


Let  H  (t ;  x)  =  P[  I 

£ 

It  can  be  seen  that  H  (t; 


S  t  j  X--x]  . 


and 


H(t) 


H(t;  x)do  (x)  , 

J0 


where  the  integrals  in  this  section  are  taken  to  be  Lebesgue*' 
St  ie  1 1 jes  i  ntegra 1 s . 

Cone!  i  t  i  on  i  n g  on  X  , 


40  h iiuidu  a  r  r  .  du 


'  tQuld 

t  HC(t) 


t“0  H  (t) 


r  r  .  du .  (X) 

J0  Hc(t)  vt  Hc(t;x) 


since  the  integral  converge  absolutely, 


.  f  Hbluxi  f  UKuM u  .  d,  (x) 
0  HC  (t)  t-i-x  Ac  (t-i-x) 


-  P  A^Xyldu  r  Hl(l;xl  d?  (*) 

(S)  Jt  AC(t)  0  HC(t) 


from  the  assumption  of  the  theorem. 


Hence , 


03  HC  (u)  du  o'  ,C°  Ac  (u)  du 
HC(t)  (£)  Jt  AC(t)  ’ 


By  letting  t~0  one  gets  the  following 


Corol  la  ry  :  For  (DMRL)/G/1  queues,  ( ( 1  MR  L)  /G/l  queues), 


-  1 


Theorem  2.2:  For  all  Gl/G/1  queues  where  A (t )  is  restricted  to  be  IFR  (DFR)  , 


■  V  .  - 


for  A  >  0  ,  and  all 


HC(t)  &  AC(t) 


t  22  0  where  f  i n  i  tc  , 


di)  1 


all  150, 


A  (t)  (t) 


OH)  I' 


>c0  1 1 c  /  \  |  ^  pc0  ,C 

H  (u) du 

t  vh  ~»J 


AlCyiiii  all 


t  va 


1  5  0. 


The  inequalities  in  (i)~(iii)  are  tight  for  the  M/G/l  queue,  and 
in  this  case  the  ratio  in  (ii)  is  constant  and  equal  to  1. 

Proof :  Proofs  are  given  for  A(t)  not  discrete. 

(i)  First,  condition  on  X  as  in  the  proof  of  the  previous  theorem. 


H(t-tA)  -  H (t)  =  r  fH(t-i-A;x)  -  H ( t ; x)  ]  d§ 
HC(t)  ^0  HC(t) 

Now  H(t+A;x)  -  H(t;x)  =  Hc(t;x)  -  HC(t+A;x) 


and  substituting  we  get 


Hj.t1: A. Hill  s 1 f  AJ^t^„-_Alt±>sl  .  d5  (x) 

HC(t)  HC(t)  J0  AC(x) 

_ ! _  r  A(t-!-A-i-x)  -  A(t-i-x)  o  AC(t+x)  d6  ^ 

HC(t)  0  AC  (t-i-x)  AC(x) 

MJJA) — from  the  I FR ,  (DFR)  assumption  on  A(t). 
(•')  AC(t) 


( i i )  Add  and  subtract  1 


from  both  sides  of  (i), 


i  _  \til3l  - 

HC(i)  <■' 


AC(tiA) 

AC(t) 


xt 
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or 


«!iii  3  >fhm  anA>0,  ts  o, 

Ac(t)  ^  Ac(t-l-A) 


which  proves  part  (ii).  Notice  that  (i)  and  (ii)  are  equivalent 
as  this  argument  is  reversible. 


(i'i) 


From  theorem  2.1 

rfO 

Hc(u)du 
Jt  == 

P“  r  (=) 

A  (u)du 

Jt 


AC(t)  (S)  AC(v) 


all  0  S  v  £  t 


from  (ii)  and  the  fact  that  IFR  (DFR)  ->  DKRL  (IMRL)  (sec  Barlow 
and  Proschan  ( 1 96f3 ) )  . 

Putting  this  in  determinant  form, 


1— '  r*1'*" 


St2£HiAxar**m:- 


Adding  the  first  column  to  the  second  gives 


<N) 


0 


which  proves  (iii)  and  completes  the  proof  of  the  theorem. 

Part  (iii)  leads  to  the  following 
Coro!  la  ry ;  For  1FR/G/1  queues  (DFR/G/1  queues) 


.(2)  J2) 


g  Va 


2v,  (?.)  2v 

h  '  '  a 


X  /  2  / 1  \  2\  'a 

2  <°a  +  )  = 


K  +  1) 


2X 


Equality  is  taken  on  everywhere  by  the  M/G/l  queue. 


3 .  Bounds  for  Queues  with  Monotone  Failure  Rate  Arrival  Distributions 
a)  The  probability  an  arrival  does  not  wait. 

Recall  from  Chapter  1,  theorem  1.1  that 

Using  the  corollary  of  theorem  2.1  above, 

(i)  For  DMRL/G/1  queues 


(1-p)  S  ao  ^  1. 


(ii)  For  IMRL/G/1  queues 


y \  SSSCgg*  TVTT-t- 


0  5  a  2=  (l-o) 
o  •  ' 


The  lower  bound  in  (i)  and  upper  bound  in  (ii)  are  taken  on  by  the 
Poisson  arrival  queue.  The  upper  bound  in  (i)  is  taken  on  by  the 
constant  arrival,  constant  service  queue, 
b)  The  mean  idle  time. 

Using  equation  11  of  Chapter  1  and  the  corollary  to  theorem  2.1  above 
gives 

(i)  For  DMRL/G/1  queues 


tfXH.)  a  »h  a  £ 


(ii)  For  IMRL/G/1  queues 


X  -  V 


The  upper  bound  in  (i)  and  lower  bound  in  (ii)  are  taken  on  by  the 
Poisson  arrival  queue.  The  lower  bound  in  (i)  is  taken  on  by  the  D/D/1 
queue . 

c )  The  mean  length  and  number  served  in  busy  period . 

Using  the  relationships  E[B]  ~  ypyy  E[l]  and  £[1^]  =  pE[B], 

(i)  For  DMRL/G/1  queues, 


-  S  E[B]  §  -7-j—T, 
k  p(l-o) 


1  -  V  5  TW) 


‘’•-a*  —  -  • 
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the  upper  bounds  being  tight  for  the  M/G/1  queue,  and  the  lower  ones 
for  the  D /D /I  queue . 

(ii)  For  1MRL/G/1  queues, 


_J _ 

pO-p) 


^  E[B], 


1 

(1-p) 


< 


E[NJ. 


These  bounds  are  tight  and  taken  on  by  the  M/G/1  queue. 

d)  The  mean  wa  i t  and  number  in  queue. 

Here  use  is  made  of  the  important  queueing  formula  E[N  ]  =  XE[W]  (see 

9 

Little,  (1961)).  With  this,  equation  (12)  of  Chapter  1  and  the  corollary 
to  theorem  2.2  above,  the  following  results  are  obtained. 

(i)  For  all  IFR/G/1  queues 


(P+Ca) 

J  ~  ~ 2T“  g  E[W]  "  J’ 


(22) 


and 


where 


(p+c?) 

Xj  _  — E[Nq]  XJ, 

2.2  2 
c  +  X  a 
.  _  _ 9L 

2X ( 1 -p)  * 


(23) 


The  first  expression  shows  that  for  this  broad  class  of  queues  the 

expected  wa i t  has  been  bounded  to  within  at  most  a  mean  inter-arrival 

2 

time  since  p  <  1  and  for  IFR  arrivals  c  ^  1  (see  Barlow  and  Proschan 

a 

(1965)).  The  second  expression  shows  that  the  expected  number  in  the 
queue  has  been  bounded  to  within  at  mos  t  1  customer. 


,t=vT^*snr- 


TTVBV^ 
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(ii)  For  all  DFR/G/1  queues 


c2  +  X2cj2 


(p  +  C  ) 
a 

2\ 


and 


s  «V  *  20-00 


2  2  2 

c  +  Xo  p  +  c 
a _ q  _ _ a 


(24) 


(25) 


where  t  solves  equation  (14).  The  upper  bound  is  taken  on  by  the  M/G/l 

queue.  It  has  not  been  shown  that  the  lower  bound  is  taken  on  by  any 

member  in  the  class. 

e)  The  variance  of  the  wait  in  queue. 

Theorem  2.2  allows  bounds  to  be  obtained  on  the  variance  of  an 

2 

equilibrium  excess  idle  distribution,  o g  ^ .  If  the  first  three  moments 
of  the  arrival  and  service  distributions  are  known  this  is  the  only 
unknown  quantity  in  equation  (5),  which  gives  the  variance  of  the  wait. 
C  lear  ly , 

(i)  For  IFR/G/1  queues, 


,(3)  J3) 


v 


0  So 


e,h  3v, 


3v_ 


and  this  with  equation  (5)  gives  upper  and  lower  bounds  on  the  variance 

2 

of  the  via  i  t  which  differ  by  no  more  than  since  if  A(t)  is  IFR 


v(3) 


(see  Barlow  and  Proschan  (1955)). 


TTOTS' 
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f )  The  variance  of  th e  output  dist  r  i  b  u  t  i  on  . 

The  only  unknown  factor  in  the  variance  of  the  output  i  r»  again 


2vh 

(0 


and  using  theorem  2.2, 
For  IFR/G/1  queues, 


2  <  u  r  a  -  2  ,  1-n 

o'  ^  Var  I  t  l=o  -i-  — r- 
9  n  g  ^2 


the  upper  bound  being  equality  for  the  M/G/l  c|ucuc,  and 
(ii)  For  DFR/G/1  queues, 


9  x2 


-e—  S  Var  [t  ]  S  2o2  +  o  j  -  21  (f  -  ~)  , 


VX  p/ 


where  l  solves  equation  (l4)„ 

9 )  The  Virtual  wa i t 

Using  the  bounds  in  (22)  and  equation  (6)  give, 
(i)  For  all  IFR/G/1  queues, 


p(ca+cf)  p(Cq+Ca)  (P'!'ca) 


2|i0-p) 


( 1 -p  ) 


2p, 


and  the  expected  virtual  wait  has  been  bounded  to  within  at  most  — . 

The  lower  bound  is  equality  for  the  M/G/l  queue.  The  upper  bound  is 
the  Cesaro  mean,  or  the  mean  wa  i  t  of  a  random  arrival  in  the  D/D/1  queue, 
(ii)  For  all  DFR/G/1  queues, 


>[ 


(Hc2/ 


i]  a  e[v] 


/  2,  2. 

P  (cci'ca> 
2n  ( 1  -p )  ' 
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4 .  A  Weaker  Assumption  on  the  Arrival  Stream 

Bounds  on  the  expected  number  in  the  queue  to  within  at  most  one 
customer  can  be  obtained  by  weakening  the  IFR  assumption  on  the  arrival 
stream.  We  first  prove  the  following 

Theorem  2.3t  For  all  Gl/G/1  queues  where  A(t)  satisfies 


I 


AC  (u)  du  =E 
t  AC(t)  ^ 


M  <  “  for  all  t  I  0, 


v 


(2) 


then 


2vh  (S) 


M. 


Equality  holds  when  A(t)  is  exponential  and  M  = 


1 


Proof :  The  proof  is  similar  to  that  of  Theorem  2.1.  Using  the  same  notation 
as  given  there, 


ot  r 
r* 


h1m.cJ£  =  r  Hfitixi  r  ^(ujdu  r,~ ^ 

t  HC(t)  ^0  HC(t)  t+x  AC (t-i-x) 


CO 


M  [  di  (x)  -  M 

J0  HC(t) 


from  the  assumption  of  the  theorem. 
Therefore  , 


all  t  1=  0. 


Integrating  over  t  gives 


GO 


M 


and  the  theorem  is  proved. 


- 


3?-. 


Using  this  result  and  (12)  in  (4)  gives, 
For  all  GI/G/1  queues  with  p  <  1  and 


f00  AC  (u)  du 
Jt  AC(t) 


X 


for  a  1  1  t 


0, 


J 


Lhul 

2X 


g  E[W]  g  J  , 


where 


EtN^'J  S  XJ, 


J  = 


2 

c 


,  -v  2  2 

+  X  o 


a _ g 

2X(!-p) 


i  Example :  Suppose  we  have  m  types  of  customers  each  requiring  m  different 

types  of  service,  and  we  have  one  facility  for  each  type  with  service 

1  2 

distribution  G.(t),  mean  — and  variance  a 
i  u, .  g. 

r  i  3 1 

Customers  arrive  in  a  renewal  process  with  inter-arrival  time  distri¬ 
bution  A(t)  which  has  DMRL  (1MRL).  The  probability  that  an  arrival  is  a 

m 

type  i  customer  is  p.  5  0,  where  £  p .  =  1 .  The  distribution  of  time 

i  =  l 

between  two  type  i  arrivals  does  not  in  general  have  DMRL  (IMRL).  For 
example  if  in  the  DMRL  case  the  original  distribution  A(t)  is  degenerate 
at  r-,  the  distribution  of  time  between  type  i  arrivals  is  a  step  function 

A. 

which  does  not  have  DMRL.  However,  we  show  next  that  the  assumption  of 
Theorem  2.3  holds  for  each  i,  when  M  ~ 

^  i 

Consider  a  renewal  process  of  a  1 1  arrivals  starting  at  t~0  with  a  type 

i  arrival.  Pick  some  arbitrary  time  tj  >  0,  and  define 

Xj  =  time  between  (j-l)-th  and  j-th  arrival  after  t^,  j-1  ,  2,  ...  . 

X  =  time  to  n-th  arrival  after  t,. 
n  >t,  1 


*5$  y 


I 


33. 


Then 


X  ,  =  T,  X. 

n,tl  j*=1  J 


E[X„  ]  =  S  E[X.3  -  E[X,]  +  (J,  2. 


n-tr  j=i  - j 


from  DMRL  (IMRL)  property.  Now  let  X  =  time  to  first  type  i  arrival  after  t. 

L  ■,  • 


Then 


E[Xt,]  |)  r  S  nPid-Pi)"-'  = 

1  v  '  n=l  i 


independent  of  when  the  last  type  i  arrived  before  t|.  Let  T.  -  time  of  last 
type  i  arrival  before  t ^ .  Then  in  particular 


|)  x; 


But 


A*:  (ix)  dp. 

EC*JTr°W 


t,  * j  (t  j ) 


and  the  result  is  at  hand. 


v 


(2) 

h. 


It  follows  immediately  from  theorem  2.3  that  —  S  M-  i  =  l  ,  2,  ...,  m 

h .  P  i 

i 


and  for  the  i-th  facility 


(1+p,) 

Ji  ■  “sp- s  E[w]  s  Ji- 


(22a) 


(1+Pi) 


Xp.J:  -  • — rC" 

r  i  i  2 


S  E[Nq3  S  Xp.J., 


(23a) 


J. 


p.c  +  (Xp.)  o  -l-  1-p. 
“  i  a  w  r  g.  r  i 

_ _ 

2Xp . ( 1 -p  . ) 


Xp 

i 


whe  re 


3*i. 

5.  Exponential  Service  Queues 

The  GI/M/1  queue  has  been  studied  at  length  by  various  authors  and  an 
important  quantity  to  be  determined  in  its  analysis  is  6,  the  unique  non¬ 
zero  root  of  the  equation 

z  =  “a  (|a  (1-z))  (26) 

-CO 

where  'a(s)  =  e  stdA(t). 

V 

It  has  been  shown  that  (for  example  see  Prabhti  (1965  o)) 

aQ  =  P[W=0]  =  (1-6), 

P[V/^x]  =  6(1  -e  ^ ^  ^ ) x )  o<x, 

P[an  arrival  finds  j  in  system  I  -  (1-6)6-*  j-0,  1,  2,  ...  . 

By  restricting  A(t)  one  can  say  something  about  the  magnitude  of  6 
in  comparison  to  p  =  In  fact  we  prove  here 

Theorem  2.4:  If  A  ( t)  is  IFR  (DFPO  and  if  6  is  the  unique  solution  of 

^  — 

{ 0  <  z  <  1  ,  z  =  a  (p  ( 1  -x) )}  for  -  =  p  <  1  ,  then  6  ^  p. 

Equality  is  taken  on  by  the  exponential  distribution. 

Proof :  The  proof  is  a  direct  application  of  the  corollary  4.9  to  Theorem 
4.8  on  page  33  of  Barlow  and  Proschari  (1965).  There  it  is  shown  that  if 
A  (t)  is  IFR  (DFR)  with  mean 


.tv 


'aESii'SSk 


I 


35. 


a'(s) 


?5  X 


(=)  X+s 


for  s  rea 1 . 


X-i-s 

X 


For  a(s)  =  the  solution  of  |o  <  z  <  1  ,  z  =  * s  z  ~  ~  "  P  • 


It 


X 


For  a(s)  Z?  ~jr  the  solution  must  lie  to  the  left  of  p  arid  for  a(s)  5  ^.ls  , 


it  must  lie  to  the  right. 

This  leads  to  the  following: 
For  the  IFR  (DFR)/M/1  queue 


1)  P[W  >  x]  =  6e"p'^1"6^X  ^  pe"^(1”p)x  ¥  x  5  0, 


2) 


E[W] 


6  -  P 

X(l-6)  (5)  X ( 1  -p ) 


3) 


E[Nql 


JL  -  „_P__ 
1-6  (S)  (1-p)’ 


Equality  holds  throughout  for  the  M/M/1  queue. 


6 .  Expected  Number  found  by  an  Arrival 

For  the  M/G/l  queue  it  can  be  shown  that  E[N  *],  the  time  average  number 
in  the  queue  is  the  same  as  the  average  number  found  by  an  arrival.  Let 
L  be  the  expected  number  in  queue  found  by  an  arrival,  and  let  L  =  E[ N  ]. 

3  t  C| 

For  the  Gl/G/l  queue  one  can  shov/  that 

Lt  -  PLa  -i-  (l-ao)XY,  (27) 

where  y  -  average  excess  service  time  of  the  customer  in  service  when  an 
arrival  comes.  The  argument  is  similar  to  that  leading  to  equation  (6)  for 
the  virtual  wait.  For  Poisson  arrivals  a  =  (1-p), 


I 


Chapter  3 

BOUNDS  FOR  BATCH  ARRIVAL,  BATCH  SERVICE, 
AND  INTERRUPTED  SERVER  QUEUES. 


37. 


1 .  I  n trod  act i on 

This  chapter  considers  some  generalizations  of  the  GI/G/1  queue  which 
preserve  some  important  basic  points  in  its  underlying  mathematical 
structure. 

V/e  first  study  a  queue  where  customers  arrive  in  batches  of  (possibly) 
random  size.  This  is  treated  by  redefining  a  service  time.  Queues  with 
services  occurring  in  batches  of  fixed  size  N  are  then  analyzed,  and 
solved  by  redefining  an  inter -arrival  time.  Finally  a  queue  in  which  the 
service  mechanism  is  subject  to  breakdown  in  a  Poisson  manner  with  a  general 
repair  time  distribution  is  studied. 

Each  of  these  three  generalizations  of  the  basic  GI/G/1  queue  pre¬ 
serves  the  following  structure.  The  (possibly)  redefined  inter-arrival 
or  service  times  both  form  independent  sequences  of  independent  and 
identically  distributed  random  variables.  Also  a  sequence  of  waiting  times 
can  be  found  which  satisfy  the  recursion  formula 


W 


n-i-1 


Max 


(0 ,  W  +S  ' ; 
'  n  n 


■T')  , 
n' 


where  the  primes  denote  a  possibly  redefined  service  or  inter-arrival  time. 


2 .  Queues  with  Batch  Arrivals 

Suppose  "arrivals"  appear  at  rate  X  in  a  renewal  process  but  that  each 
"arrival"  is  now  a  batch  of  customers  of  random  size  T|.  Let  P(T)-n)  -  rr^  , 


38. 


n=0,  1,  2,  ...  .  Customers  in  each  batch  are  assumed  to  be  numbered  in 


some  way  to  denote  their  order  of  service.  The  FIFO  order  of  serving 


ba  tches  is  maintained.  Note  that  we  allow  a  zero  batch  size,  which  is  of 


significant  importance  as  will  be  shown  in  an  example. 


Now  let  S.  be  the  service  time  of  the  k-th  customer  in  the  n-th 
k,n 

batch  of  a  stationary  queue,  and  S,  ~G(t)  for  all  n  and  k=  1,  with 

K  j  n 

Sq  n  =  0.  Let  V/^  n  be  the  waiting  time  (in  queue)  of  the  comp  1  e  t  e  n-th 

batch,  and  define  S'  =  Sn  +  S,  +  . . .  +  S»  ,  S*  ~  G*  (t)  .  Then  the 

’  n  0,n  l,n  T]  ,n  n  7 

sequence  satisfies  the  relationship, 


W.  .  .  =  Max  [0,  W.  +S" -T  1 , 
b  ,n+ 1  b  ,n  n  nJ 


where  T  is  the  time  between  the  arrival  of  the  n-th  and  (n-M)-th  batches, 
n 

T  ,~A(t)  for  all  n  as  before, 
n  ' 

Noting  that  {S^}  and  fT^}  form  independent  sequences  of  independent 
random  variables,  when  \v  <  p  the  theorems  in  Chapter  1  are  valid. 

The  expected  wa  i  t  of  the  first  customer  in  each  batch  is  thus  given  by 


2.  2 
a  +a  , 
a  q* 

v 


v 


2<r 


TTs 


,  1(1  - 
2h  n 


TT 


V 


(2) 


■)  - 


2v. 


(30) 


V> 


The  variance  a  ,  can  be  found  easily  by  conditioning  giving 
fp' 


a 


2 

g* 


2 

v  a 
tt  g 


+ 


2 

a 

IT 

2* 

V> 


(31) 


The  average  waiting  time  of  an  arbitrary  customer  is  found  as  follows. 
The  expected  total  additional  wa i t  of  a  1 1  customers  in  an  average  batch  is 
found,  and  this  is  divided  by  the  average  number  of  customers  per  batch. 


*UtS 
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Let  X  be  the  total  additional  waiting  time  of  all  customers  in  some  batch 
n 

n.  Then  X  -0  if  T)=0 ,  1,  and  for  T|  S  2, 
n 


X  =  S.  + 
n  1  ,n 


(S.  +S-  )  + 

'  1  ,n  2  ,n' 


+  (S  .  +S  0  + 

'  1  ,n  2 ,n 


+s  (H-l)  ,nh 


Conditioning  on  'H  it  follows  that 


E[X  I'm  =  1 ^  for  all  T)  S  0, 

n 1  "  <-p, 


and  hence 


E[Xnl 


(2) 

V  '  '  -V 
rr  tt 


2^, 


The  expected  wait  of  any  unspecified  customer  E[V/1  is  then  given  in  terms 
of  the  wait  of  the  first  one  in  each  batch  by 

v(2) 

E[vn  .  E[w  1  ♦[£--£]  -■ 

TT 

Using  this  with  (31)  and  (30)  an  expression  is  obtained  for  the  expected 
wait  which  is  bounded  as  follows  under  the  assumption  of  1FR  arrivals 
(between  batches): 


(o+c^) 

VJ2  --IT-SEOIS  J,+J2 


v(p+c) 

Xv  J,  +  Xv  J_  -  ^  e[N  1  g  Xv  J  +  Xv  J 

tt  1  tt  2  2  q  J  tt2  tt 


(32) 

(33) 


whe  re 


and  J 


,r  2  2  2  2  2-, 
AUl  cr  -In  v  O'  -la 
1  a  _  TT  q  Tf 

2h2(I-p) 


Xv 


TT 


V> 


-V 


(2) 


TT 

,2v 


TT 
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Hence,  the  average  number  of  customers  in  the  queue  is  bounded  to  within 
at  most  an  average  batch  size.  For  the  compound  Poisson  input  the  lower 


For  queues  where  f  - — (iiL-L- ^  ~  for  all  t 

Jt  Ac(t)  X 


bound  is  tight 


5  0,  we  can  replace 


p-i-c 


w.  ,  (o  +  l) 

lT  by  2X  • 


It  is  of  interest  to  note  that  the  bounds  still  hold  even  though  the 
distribution  of  time  between  the  arrival  of  cus  tomers  is  not  IFR.  For  a 
random  variable  to  have  an  ! FR.  distribution,  the  distribution  can  have  at 
most  one  jump,  that  being  at  the  right  hand  end  of  its  interval  of  support 
(for  details  see  Barlow  and  Proschan  (1965)).  It  is  obvious  that  if  there 
is  a  positive  probability  of  having  two  or  more  customers  in  a  single  batch 
then  there  is  a  non-zero  probability  of  a  zero  inter-arrival  time  between 
customers.  In  fact  let  T"  be  the  time  between  customer  arrivals,  T  ~A  (t)  , 
then  considering  only  batches  of  non-zero  size  we  find  that 


A“(0")  =  0  and  A”(0'!)  = 


V  -  (  1  -TT  ) 

TT  _ cr_ 

v 

TT 


These  inter-arrival  times  are  neither  independent  nor  IFR. 

As  a  special  case  of  batch  arrivals  consider  the  example  in  section  h 
of  Chapter  2,  where  a  customer  from  the  original  arrival  stream  stops  at 
facility  i  with  probability  p..  Since  in  our  development  we  allow  batches 
of  zero  size  we  can  interpret  this  example  as  batch  arrivals,  where  the 
batch  size  is  1  with  probability  p.  and  0  with  probability  ( 1 -p . ) .  Hence, 


v  =  p.  and  o  =  p.(l-p.).  From  (32)  and  (33),  for  the  case  where  A(t) 
rr  i  nil 

is  assumed  to  have  bounded  mean  residual  life  (as  in  section  h,  Chapter  2), 


rP  ;  + 

j.  -[irJsEM «  j,. 


(32a) 


Xp .  J .  - 

'  i  i 


P:  (P  ;+l) 


-  E[NJ  *  Xp:J:, 


(33a) 


whe  re 


X  = 


p.c^  -i-  (>p.)^a^  -i-  p?(l-p.) 
“  i  a  v  i _ q  H  i  v  r  i 

2Xp .  ( 1  -p  . ) 


Comparing  thes  .•  with  (22a)  and  (23a)  we  see  that  this  approach  gives  a 
better  bound,  as  the  expected  number  in  queue  is  bounded  to  within  at  most 
p.  customers.  It  is  easy  to  show  the  lower  bounds  are  the  same  in  both 
cases  and  are  taken  on  by  Poisson  arrival  queues. 

It  is  also  important  to  note  that  the  batch  arrival  approach  to  this 
example  gives  even  tighter  bounds  if  A(t)  is  assumed  to  be  IFR,  whereas  the 
approach  in  Chapter  2  breaks  dov/n  in  this  case,  since  as  shown  above  the 
distribution  of  time  between  customer  arrivals  is  not  IFR.  The  above 
approach  gives  the  expected  number  in  queue  to  v/i  thin 


Pi  (p  +o 


if  A ( t )  is  IFR  . 


All  the  DFR  bounds  of  Chapter  2  hold  when  the  correct  means  and 
variances  are  defined,  since  the  idle  time  is  still  some  excess  of  an  inter¬ 


batch  time. 


3 •  Queues  With  Batch  Services 

In  this  section  we  assume  a  renewal  input  into  a  single  channel  queue 
v/i  th  independent  service  times,  but  where  service  takes  place  only  in 


batches  of  a  fixed  size  of  N  customers  ervice  and  i n ter -a rr i va 1  times 
are  again  mutually  independent. 

Let  ^  be  the  wait  in  queue  of  the  last  (N-th)  person  who  arrived 
and  formed  the  n-th  batch  of  N  customers.  By  redefining  an  inter-arrival 
t i me  to  be 


A  ^  ^ 


where  T.  is  the  i  n  ter-a  r  r  i  va  1  time  between  the  j-th  and  (j-i-l)-th  members 

J.n 

of  the  n-th  batch,  we  have 


Wm  x i  =  Max  [0,  WM  +  S  -T“], 
N  ,n+ 1  L  ’  N,n  n  nJ 


where  S  is  the  service  time  of  the  n-th  batch, 
n 

The  new  inter-arrival  time  distribution  is  the  N-th  fold  convolution 
of  the  original  A(t).  The  IFR  property  is  closed  under  convolutions  and  so 
with  this  assumption  on  A(t)  the  arguments  of  Chapter  2  can  be  applied. 
However,  the  DFR  arguments  fail  since  convolutions  of  DFR  distributions 
are  not  necessarily  DFR  (for  details  see  Barlow  and  Proschan,  (1965)). 

For  the  distribution  A*  (t)  we  see  that 


v 


N 

X  ’ 


a  ,  =  No 
a*  a 


and 


For  any  unspecified  customer 


E[W]  =  E[W 


+ 


1  _0iiD_ 

X  2  ’ 


N 


and  hence  for  IFR/G/1  queues  with  service  in  batches  of  fixed  size  N 


whe  re 


(Np-i-c^) 

J  "  —2)T~  "  E[W^  2  J’ 


2  2  2 
Nc  -I-  X  a 

a _ g. 

2  NX  ( 1  -p )  * 


(34) 

(35) 


Thus  the  expected  number  in  the  queue  is  bounded  to  within  at  most 
cus  tomers . 

It  is  well  known  that  in  generalizing  the  M/G/l  queue  to  batch  services 

analytic  problems  arise  and  essentially  the  problem  is  the  same  as  analyzing 

the  E^/G/l  queue.  In  order  to  calculate  E[\7l  or  E[N  ]  for  this  case  it  is 

necessary  to  determine  the  roots  of  the  equation  (s-i-X)^  =  X^g(-s)  for 

00 

Re(s)  l-  0,  where  g(s)  =  e  J'dG(t)  (for  details  see  for  example  Prabhu 


0“ 

1965  (a)).  However,  all  this  can  be  avoided  if  the  approximations  in  (34) 
and  (35)  give  sufficient  accuracy. 


4 .  Queues  wi th  Serve r  Breakdown 

We  now  consider  a  GI/G/1  queue  whe  re  the  service  mechanism  can  break 
down.  V/ hen  this  occurs  the  normal  queueing  process  is  interrupted  until 
the  server  is  repaired,  and  it  is  assumed  that  these  interruptions  occur 
only  in  busy  periods . 

It  will  be  assumed  that 

a)  Times  between  successive  breakdowns  of  the  server  are  exponentially 

distributed  and  independent  with  mean  ^ — . 

^b 


Vi. 

b)  Repair  times  are  independent  drawings  {R^}  f r o:n  a  distribution  R(t). 

c)  The  customer  in  service  when  an  interruption  occurs  resumes  his  service 
where  he  left  off.  This  is  called  the  Resume  Rule. 

To  analyze  such  a  queue  we  use  a  concept  introduced  by  Gaver  (1962) 
and  dealt  with  in  detail  in  his  papers.  This  is  the  "completion  time"  of  a 
customer,  for  any  customer  n.  It  is  the  time  the  n-th  customer  spends 
actually  occupying  the  service  facility  (his  service  time  and  all  the  repair 
time  for  breakdowns  during  his  service).  The  assumption  of  exponential 
i nter -breakdown  times  leads  to  {Cn}  being  a  sequence  independent  and  iden¬ 
tically  distributed  random  variables.  The  sequence  of  waiting  times  (in 
queue)  satisfy 


W 


n-l- 1 


Max  [0,  W  -l-C  -T  ]. 

n  n  n 


(36) 


For  the  IFR  or  DFR  assumption  on  the  arrival  stream  all  the  bounds 
of  Chapter  2  hold  in  terms  of  the  mean  and  variance  of  the  inter-arrival 
times  and  the  mean  and  variance  of  a  completion  time  C.  These  last  two 
quantities  are  easily  shown  to  be 


E[C]  =  -HI  -tv  A,  ) 
J  p  r  b 


V[C  ]  =  (1  +  v  A  )V  -t  —(a2  tv2) 

L  J  '  r  b'  g  p  v  r  r 


For  a  stationary  queueing  distribution  to  exist  we  must  have  y  >  E[C]. 

For  any  other  distribution  of  time  between  breakdowns  the  completion 
times  will  not  in  general  be  independent  or  identically  distributed  and  the 


method  fa i  Is. 


We  can  replace  assumption  c)  by  c^), 

C|)  The  customer  in  service  when  an  interruption  occurs  starts  his 

service  over  again,  but  with  the  same  service  time.  This  is  colled 
the  Repea  t  I  den  t i ca 1  Rule. 


or  by  c^) 

c^)  The  customer  in  service  when  ar.  interruption  occurs  starts  over 

again  with  a  new  independent  service  time.  This  is  called  the  Repea i 
D  i  fferen  t  Rule  . 

For  a  detailed  discussion  of  these  see  Gaver  (1962). 

In  both  the  above  cases  the  completion  times  form  sequences  of  inde¬ 
pendent  and  identically  distributed  random  variables  which  are  independent 
of  the  arrival  stream.  Hence,  all  results  for  the  Gl/G/1  queue  follow 
since  consecutive  v/aiting  times  satisfy  the  fundamental  equation  (36). 

Under  the  IFR  or  DFR  arrival  assumption  the  relevant  bounds  are  thus  ob¬ 


tained.  Also  under  the  assumption  that 


^  AC(u)du 


the  relevant 


t  AC(t)  X 

bounds  follow.  All  give  bounds  on  the  expected  number  in  queue  to  within 
at  most  1  customer. 
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Chapter  4 

QUEUES  WITH  ADDED  DELAY  FOR  THE  FIRST  CUSTOMER 
IN  A  BUSY  PERIOD 

1  .  I n  troduct i on 

V/e  now  consider  a  single  channel  queue  with  renewal  input  at  rate  >  , 
independent  service  at  rate  p, ,  but  where  the  first  customer  in  every  busy 
period  suffers  a  random  delay  R  before  his  service  commences.  This  random 
variable  R  can  be  dependent  on  the  preceding  inter-arrival  time  and  may 
affect  the  independence  of  the  wa i t  in  queue  of  some  customer  n  and  the 
inter-arrival  time  of  the  n-th  and  the  (ni-l)-th  customers.  It  is  assumed 
that  R  is  independent  of  the  service  times  of  all  customers  in  the  busy 
period  generated  by  R. 

In  this  chapter  the  random  variable  I  will  refer  as  before  to  the  time 
from  the  end  of  a  busy  period  until  the  next  customer  arrives.  Hence,  it 
is  still  some  excess  of  an  in  ter -a r r i va 1  time.  The  idle  time  with  respect 
to  the  server  will  be  I -i-R . 

A  general  expression  is  found  for  the  expected  waiting  time  in  queue 
and  is  seen  to  be  a  generalization  of  equation  (4)  in  Chapter  1.  As  an 
example  v/e  study  the  single  server  queue  where  the  server  waits  until  m 
people  are  present  before  starting  on  the  first  service  in  a  busy  period. 

He  then  continues  to  serve  one  customer  at  a  time  until  the  system  is  empty. 
V/c  shall  call  this  the  G I  /G  / 1  queue  (see  Heyman,  I96G)  .  Vi  hen  m~l  the  sub¬ 
script  is  dropped.  For  IFR  arrivals  bounds  ore  obtained  which  give  the 
expected  number  in  queue  or  system  to  within  ~  customer.  For  arrivals  with 

mean  residual  life  bounded  above  by  ~  the  expected  number  is  bounded  to 
2 

-nrt  1  -c^-j 
w  i  thin  - — 2n  ~C~J  • 


2.  General  Results 


By  equating  input  with  output  we  shall  n ow  find  an  expression  for  the 
expected  waiting  time  in  queue. 

Let  R  be  the  delay  before  service  commences  for  the  first  customer  in 

a  busy  period  in  a  stationary  queue.  This  may  be  dependent  on  the  previous 

inter-arrival  time  and  hence  on  the  idle  time.  It  may  also  affect  the 

independence  of  V/  and  as  is  seen  in  the  example  in  section  3  (in  this 

example  the  wait  of  the  first  customer  in  a  busy  period  is  in  fact  equal 

to  the  next  (m-1)  inter-arrival  times).  It  is  assumed  that  it  has  no 

effect  on  the  independence  of  17 n  and  S  .  \7riting  the  inter-output  times 

as  t  =  S  , ,  +  X  ,  when  customer  n  leaves  the  system  busy,  X  -0,  and  when 
n  n-l-1  n  ’  ’  n 

he  leaves  it  empty,  -  |-i-R,  independent  of  S  p 

We  now  show  that 

ET  U21  E[R?'l  -  Ef  l?'l  C°V  (V/n’Tn 

“  -2E[UT  2(eIr]  '+  E[TiT  ECU] 


Note  that  this  is  a  generalization  of  equation  (4)  in  Chapter  1. 

To  prove  (37)  we  proceed  as  follows: 

Equating  input  and  output  times 


•Eanfci®.':. 


I 


1 


where  a) 

b) 

c) 

Squaring  both  sides  of  (38),  taking  expectations  and  assuming  stationarily 
in  the  queue,  we  get 

-2E[U]E[W]  =  ECU2]  +  E[Xn(Xnn-2(Dn-Tn)]  -  2Cov  (Wp  ,Tn)  .  (39) 

Using  a)  and  b)  above  we  see  that 

E[Xn(Xn+2(Dn-Tn)]  =  E[(R+l)(R-l)]p  (40) 

where  p  =  P[  (D  -T  )  <  0], 

n  n  J 

But  from  c) 


48, 


D  -T  £  0  =>  X  =  0, 
n  n  n 


D  -T  <  0  =>  (D  -T  )  =  -I  and  X  =  HR, 
n  n  '  n  n'  n 


E[Xn]  =  -E[U] 


(from  (38)) 


PE[R+I]  -  E[Xp]  -  -E[U].  (41) 

Substituting  (40)  and  (41)  in  (39)  gives  (37). 

Note  that  no  assumption  had  to  be  made  concerning  the  independence 
of  R  and  1,  and  if  R  is  independent  of  the  arrival  stream,  Cov  (W  ,Tn)  -  0. 

3.  The  Gl  /G/l  Q.ueuc 
- n1 - - - 

Suppose  the  server  does  rot  start  service  until  rn  customers  arc  present. 
In  this  case  the  first  customer  served  in  a  busy  period  waits  for  the  next 
(m  —  1 )  customers  to  arrive  before  commencing  service.  Numbering  the  first 
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customer  in  the  busy  period  I  for  convenience,  R  =  +  •••  H‘  ^ni  |>  So 

that 


E[R]  =  EfR2]  =  [c2  +  m-l].  (<I2) 

X 

For  this  case  it  is  obvious  that  for  the  first  (m-l)  customers  in 
each  busy  period  the  waiting  time  of  a  customer  depends  on  the  following 
inter-arrival  time,  and  we  must  now  calculate  the  non-zero  covariance  term. 

The  calculation  of  the  covariance  is  achieved  by  conditioning  on  the 
position  of  service  in  a  busy  period.  Obviously,  a  busy  period  serves  at 
least  m  customers.  We  now  use  the  covariance  relationship  for  any  three 
random  variables  X,  Y,  Z, 

Cov  (X  ,  Y )  =  E[Cov(X  ,Y|Z)]  +  Cov(E[x|z],  E[y|z]).  (43) 

We  interpret  X  to  be  a  waiting  time,  Y  the  following  inter-arrival 
time  and  Z  the  position  of  the  customer's  service  in  a  busy  period.  Note 
that  E[Y|Z'J  -  E[T.|z~i]  =  ^-,  a  constant.  Hence  in  our  case  (43)  reduces  to 

Cov(Wn,Tn)  =  E[Cov  (W  ,t|  Z)’J  .  (44) 

If  i  indicates  the  position  of  service  in  a  busy  period,  for  i=l ,  2,  . . . ,m- 

W.T.  =  (T.+T.l1+  ...  +T  .+S.+  ...  +S.  , )T. , 
ii  i  i -I- 1  m- 1  1  i  -  r  i 


and  so 


50, 


Cov(W,T|z=i)  a! 


i  =  l  ,  2  ,  ....  m-1 


=  0 


i  s=  m. 


Averaging  over  busy  periods, 


p[z-i] 


E[Nb'J 


and 


P[Z  5  ml 


i  =  l ,  2 ,  . . . ,  m- 


E[Nj-(n,-l) 


E[Mb'j 


Using  these  results  in  (44)  gives 


(m- 1 )c‘ 

Cov  (Wn  ,Tn)  - 


But 


E[Nb]  “  ’  50  finolly 


Cov  (W  ,T  )  = 


(m-1) (l~o)o\ 


n  ’  n  XE[  I’J  +  m-1  * 


(45) 


Using  (45)  and  (42)  in  (37)  gives  for  the  Gl^/G/1  queue  with  p  < 


E[WJ 


2,2  2 

c  +X  a  , 

_ _a  ,  j_ 

2X ( 1  — p )  2X 


(l-o)  + 


(m-l)[m-l-Cg1  -  X2E[  l  21 
2X"(XE[  I  ]  -I-  m-1) 


(46) 


This  expression  depends  only  on  the  first  two  moments  of  the  arrival, 
service,  and  idle  distributions,  and  on  m. 

In  general  the  distribution  of  I  will  depend  on  m  (not  the  case  for 
Poisson  arrivals).  For  any  stable  queue  (o  <  I)  I  cannot  be  identically 
zero  since  busy  periods  end  with  probability  I.  However,  by  using  the 


2-, 


apparently  crude  de  vice  of  setting  E[l]  -  E[  I  ’J  ---•  0  in  (46)  we  obtai 


i  n 


•"•EOTT'JSSTvj wwmggirar r*—  -*vwr  i 
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an  upper  bound  for  all  GI^/G/1  queues  with  p  < 


2.2  2  /  ,  2\ 
c  +X  o  /  ,\  (o-i-c  ) 

Ffi,n  -=  _£ _ a  smzLL  _  _ _ JL 

tLWJ  ::r  2X  ( 1  -o )  '  2X  2\  2X 


(47) 


We  now  show  that  this  is  actually  a  supremuin  for  the  class  of  queues 
considered  by  exhibiting  a  family  of  queues  in  the  class  with  E[l]  and 

r  2-, 

E [.  I  J  arb  i  trar  i  ly  sma  1  1  . 

Consider  a  queue  with  deterministic  arrivals  at  (normalized)  rate  1, 

1  -c 

and  deterministic  service  with  service  time  for  each  customer  - - ,  where 

m 

m  i?  2  and  c  >  0  is  arbitrary.  Then  at  stationarity  it  is  easy  to  show  that 

I  1  e 

I  =e  .  For  example  take  rn~3  ,  then  The  true  E[W]  is  easily  cal- 

|ii  J  J 

4  e  4  e 

culated  to  be  ~  -  The  above  bound  gives  a  value  of  j  +  By  letting 

e~0  in  this  example  the  bound  gives  —  whereas  the  true  E [W]  -  1,  since  in 


this  case  I  = 


3 


For  IFR  arrivals  we  now  find  a  lower  bound  on  E[Wl .  For  any  value  of 
m  it  is  easy  to  see  that  theorems  2.1  and  2.2  still  hold  for  the  random 
variable  1.  In  the  proofs  of  those  theorems  m  would  affect  the  distribution 
$  ,  (of  time  from  the  last  customer  entering  a  BP  until  the  end  of  the  BP). 
However,  the  results  were  independent  of  the  form  of  ?•  ,  and  we  can  say 
immediately  that  for  1 FR  (DFR)  arrivals  and  any  m=  I, 


E[|Pj  E[Tn]  all  n 


,  Ef  |21  ^  Efl^l 

3n  2E‘[l]  (r-i)  2E[T‘J* 


(48) 


V/e  now  prove  a  useful  leimva. 


Lemma  4.1:  If  (48)  holds  and  a  and  b  are  any  non-negative  numbers, 


3?.. 


a  -  Ef  1 21  r  a_~  _E[ T?H 
b  -i-  E[fJ  (:  )  b  -i-  Efrj  • 


P roof :  For  a  ??  0,  using  (48)  we  have 


a  S  a 

ELI]  (•=)  E[T] 


-Ef I2l  ^  -E[T21 

E[I'J  (w)  E[T]  • 


Adding 


A 1  so 


these  and  clearing  the  fraction  gives 
E[T](a-E[|2])  ^  E[l](a-E[T2]). 
b(a-E[|2])  b(a-E[T2])  for  a  ,b  ^  0. 


Adding  these  gives  the  desired  result. 

2 

V/e  apply  the  lemma  to  (46)  with  a  -  (m-l-c  )(m-1)  for  IFR  arrivals 
2  2 

since  in  this  case  c  ~  l.  For  DPR  arrivals  c  can  be  arbitrarily  large, 

a  a 

2 

so  to  insure  validity  of  the  conditions  of  the  lemma  we  let  a  =  (m-1)  in 
this  case. 

Apply  the  lemma  to  (46)  and  using  (47)  we  have  shown  for  all  IFR^/G/1 
que  ues  w i th  p  <  1 , 


J  5?  E[W]  <  J  -I-  (~)  , 


XJ  S  E[N  ]  <  XJ  ■!'  j, 


(**9) 


whe  re 


J  - 


2.2  2  ,  ,  2. 

c  +X  a  (m-l-p-c  J 

a  _a _  a 


2X ( l -p ) 


2X 


$ 
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The  lower  bound  is  taken  on  by  Poisson  arrivals,  and  as  shown  previously 
the  upper  bound  is  a  supremum  for  at  least  one  family  of  deterministic 
queues . 

For  the  weaker  assumption  of  bounded  mean  residual  life,  i.e., 


'OT  AC  (u)  du 
t  AC(t) 


is  easy  to  show  that  the  bounds  are 


J 


XJ 


(1-c2) 
_ _ a ' 

2Xm 


) 


2m 


E[vn  <  j  -i- 

E[Nql  <  XJ  + 


(49a) 


For  DFR  /G/l  queues  with  p  <  1, 
m 


c2+X2o2 

^ a  iro^r 


2X 


1 

2Xm‘ 


Equality  occurs  for  Poisson  arrivals  only  when  m=l. 


f 
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Chapter  5 

QUEUES  WITH  2  PRIORITY  CLASSES 


1  .  Arrivals  in  2  Independent  Renewal  Processes 

The  following  queueing  situation  is  assumed  in  this  section. 


1.  High  priority  customers  (called  type  1)  arrive  in  a  Poisson  manner  at 
rate  Their  service  times  are  drawn  from  some  general  service 

distribution  G,(t),  with  mean  ~ — . 

I  | 


2. 


Low  priority  customers  (called  type  2)  arrive  os  a  renewal  process  with 

distribution  of  inter-arrival  times  A_  (t)  ,  with  mean  r— \  Their  service 

A.  A 


times  are  drawn  from  a  common  distribution  G2 ( t ) ,  with  mean 


1 


1^2 


3.  The  preemptive  resume  rule  applies  to  all  type  2  customers,  (see 
Chapter  3,  section  4,  and  Gaver  (1962). 

For  convenience  in  this  section  we  shall  call  the  above  3  assumptions^. 

A  busy  period  can  start  with  the  service  of  either  a  type  1  or  a  type 
2  customer.  It  is  important  to  note  that  those  starting  with  a  type  2 
customer  have  the  same  structure  as  the  busy  periods  (BP's)  for  the  inter¬ 
rupted  server  case  treated  in  section  4  of  Chapter  3>  with  =  Xj  and 
R  (t)  the  distribution  of  a  BP  in  an  M/G/l  queue  with  the  service  distribution 
G 1  (t)  . 

For  the  study  of  this  queueing  system  we  shall  use  the  following 
notation.  A  BP  started  with  a  type  i  customer  will  be  called  a  type  i 
BP.  The  phrase  "a  type  1  service  period"  refers  to  the  time  the  server 
is  busy  serving  consecutive  type  1  customers  (i.e.,  a  normal  BP  in  the 
imbedded  M/G/l  queue  formed  by  the  type  1  customers). 


^  "V.  *"J" 


—  -  —  — -  .  ^ 
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B  -  length  of  an  a  rb  i  tra  ry  BP. 

1  -  length  of  an  arbi trary  idle  period. 

B^  =  length  of  a  type  i  BP,  i-1,  2. 

=  total  number  of  customers  served  in  a  type  i  BP,  4-1,  2. 

N j =  number  of  type  j  customers  served  in  a  type  i  BP,  i , j  —  1  ,  2. 

C  =  completion  time  of  a  type  2  customer  (see  Chapter  3,  section  4). 

h  _  ^2  i  _  r  xi  1 1_  r  x2  i  i 

Pl~[iv  P2  _  p»2  ’  K  ~  '  ^2  ’  p  Lx^J  n,  LXj+X  J  p2’ 


p 


X 

4 


and  hence  p  =  P  | '‘'P 2  * 


Arbitrary  Idle  and  Busy  Periods 

Using  the  above  definitions  and  relating  input  to  output,  we  have  for 

p  <  1  , 


E[B]  = 


E[l]. 


By  assuming  A^(t)  has  IFR  or  DFR  we  now  obtain  bounds  on  E[l]  and 
hence  on  E[B] . 

Theorem  4.1:  Under  J{  ,  if  A^(t)  has  IFR  (DFR)  and  p  <  1,  then 


E[l] 


N  ~  X, 


where  a^  (s) 


e“StdA,(t)  . 

0~ 


A2(t) 


Equality  is  taken  on  when 


;n  ninny 


P  roof :  Let  I  ~  H(t),  Y  be  the  time  from  the  end  of  the  last  BP  to  the  next 
arrival  of  type  1,  and  Z  be  the  time  from  the  end  of  the  last  BP  to 
the  next  arrival  of  type  2.  Then  Y  ~  Exponent ia  1  with  mean 

X1 

Let  Z  ~  F(t) . 

Under  the  I FR  (DPR)  assumption  we  know  from  Theorem  2.2  part  (ii) 
that  FC(t)  A'(t). 

Now  l  --Min  (Y,Z)  and  hence 


-X .  t  <  -X ,  t 

HC(t)  =  e  FC(t)  e  A^(t)  . 


Integrating  both  sides  gives  the  desired  result. 

Let  =  P(an  arbitrary  arrival  finds  the  system  empty).  Note  that 
a0  's  n°L  the  probability  an  arbitrary  customer  does  not  have  to  wait  in 
queue.  Now  using  theorem  1.1  with  theorem  4.1  and  a^  g  1  gives: 

For  all  2  priority  queues  where  J\  holds  and 
i)  A2  (t)  if  I  FR  , 


I  ■£,(*,) 


t  p  ^  o  i  )-ir  - 

-  §  E[B]  g  [ - -  r~  . 

ix  L  X,  JLl-pJ 


The  upper  bounds  are  tight  for  Poisson  type  2  arrivals. 


i i )  A„  (t)  is  DFR  , 


t-a2(X,) 


E[l], 


t 


\ 


i — 

■  i 

! _ 

Li-pJ 

Lx,  J 
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The  bound  is  tight  for  Poisson  type  2  arrivals. 

Bounds  for  the  expected  number  served  in  an  arbitrary  BP  follow 
immediately. 

Type  2  Busy  Periods 

As  noted  above  these  BP's  have  the  same  structure  as  those  in  the 
"server-breakdown"  case  dealt  with  in  section  4  of  Chapter  3.  From  the 
results  there  it  is  easily  seen  that 


E[C]  =  -~Tr1-~T. 

P2(1“P1) 

How  if  ls  is  the  idle  time  in  the  server  breakdown  case  it  is  easy  to  show 
that 


ecb(2)]  =  nrpTE[g. 


(2) 

and  we  immediately  obtain  the  following  bounds  on  E[B'  ] 
i)  For  A^(t)  IFR  and  under^  : 

(I'P2)  '  SE[B<2h  S  ’ 


(1-P)  P2 


(1  -p )  * 


(50) 


The  upper  bound  is  taken  on  by  Poisson  type  2  arrivals, 
ii)  For  A2(t)  DFR  and  under/^  : 


— ! - c  pr n  (2)n 

H2 ( 1 -0 )  -  ELB  ] 
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wi 


i  th  equality  for  ^(t)  =  !  -e 


-X2t 


Since  E 


(2) 


3^2  (1_P])  bounds  on  E^N^ 


fol low 


immediately. 


Each  interruption  of  a  type  2  service  results  in  a  type  1  service 
period.  In  each  of  these  type  1  service  periods  the  expected  number  served 

1 


is 


( 1  “P  ! ) 


Therefore,  the  expected  number  of  type  1  customers  served  in 


each  type  2  completion  time  is 


1 


P'2  ( 1  “P  ] ) 


Hence  E 


r«!a]  ■ 


w»i 


^2  ^  "P  1  ^ 


~,  and 


final  ly  we  have 


E[N(2)j  =  (n2(l-p,)  +  X,)E[B(2)]. 

(2)  (2). 

Bounds  now  follow  on  E[N'  ']  from  the  above  inequalities  on  E[  B v  '] . 

Waiting  T imes 

The  low  priority  customers  have  no  effect  on  the  high  priority  type 
when  the  pre-emptive  rule  applies.  Hence,  the  waiting  time  of  the  type  1 
is  given  by  the  Pol laczek-Kh i ntch  ine  formula.  We  can  view  type  2  customers 
in  a  similar  way  to  Chapter  4.  The  type  2  customers  form  a  general  arrival, 
general  service  queue  where  service  times  are  redefined  to  be  completion 
times,  and  the  first  customer  in  a  busy  period  has  some  random  delay  R, 
which  is  either  zero  or  the  remainder  of  a  busy  period  of  type  1  customers. 
Hence,  formula  (37)  with  Cov  (Wn  ,T  )  =  0  gives  the  expected  v/a  i  t  of  a  type 
2  customer.  Using  this  and  the  lemma  in  Chapter  4  we  get  for  the  type  two 
customers  if  ( t )  is  IFR, 
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e[(t(2)-c)2JL  +  -  e[t(2)J,  a  E[w1 

2E[T(2)-C]  2[E[R3+E[T  2  3] 


Er(T(2)-C)21  EfR21 

2E[T^-C]  2E^ 


(50 


2 

However,  to  bound  E[R]  and  E[R  ]  seems  to  be  very  difficult.  For 
Poisson  type  2  arrivals  it  can  be  shown  that,  if  R  ~R(t),  and  has  density 
r  (t)  ,  then 


RC(t)  «  B(j(t)m(X2)  -  “  r  (t)  , 

where  B^t)  is  the  distribution  of  a  busy  period  of  type  1  customers, 

00 

^(s)  |  e"StdM(t)  and  M(t)  is  the  renewal  function  for  the  modified 

J0 

renewal  process  {X.},  where  ~  exp  (X^)  and  for  i  >  1,  X.  ~  Busy  cycle  of 
type  1  customers.  Using  this  we  get 

v<*> 

ErR2i-Er  i2i  _  Jii_  j_ 

2[E[R]+E[I]]  2vb  ~  X2* 

Substituting  this  in  (37)  gives  an  exact  expression  for  E[V/]  which  checks 
with  results  inV/ei  Chang  (1965). 


2 .  Arrivals  in  a  Single  Renewal  Process 

In  this  section  we  assume  that  inter-arrival  times  form  a  single 
renewal  process.  V/ith  probability  p  the  arrival  is  a  type  1  (high  priority), 
and  with  probability  (1-p)  it  is  a  type  2.  Assuming  the  pre-emptive  resume 
rule  applies  the  completion  times  of  type  2  customers  are  no  longer  inde¬ 
pendent  and  identically  distributed  (unless  the  arrival  stream  is  Poisson). 
However,  idle  periods  (when  the  facility  is  empty)  are  excess  inter-arrival 
times  and  under  the  IFR,  DFR  or  bounded  mean  residual  life  assumption  these 
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can  be  bounded.  As  in  sec  Li  on  1  we  then  get  bounds  on  the  expected  length 
of  an  arbitrary  busy  period. 

V/a  i  t  ing  Times 


Under  the  pre-emptive  rule  the  average  wait  of  a  high  priority  customer 

can  be  bounded  under  the  various  assumptions  on  the  arrival  stream.  For 

example,  (32a)  applies  directly  when  the  inter-arrival  times  have  bounded 

mean  residual  life.  This  surprisingly  gives  a  better  bound  than  the 

approach  in  Chapter  2  (see  22a).  in  the  non-preempt i ve  case  the  type  1 

customer  w  ho  starts  a  type  1  service  period  must  v/a  i  t  until  any  type  2  in 

service  leaves  the  system.  Hence,  (31)  holds  where  in  this  case  R  is  the 

remaining  service  time  of  a  type  2  if  one  is  present,  and  R=0  if  a  type  1 

arrival  finds  the  system  empty.  If  the  type  2  customers  have  a  scrv  ice 

2. 

distribution  which  has  bounded  mean  residual  life,  then  0  5  f E^rI"  "  p~ 

and  using  this  in  (51)  gives  the  expected  v/a  i  t  of  a  type  1  customer  to  within 

r— +  — ,  and  hence  the  expected  number  of  type  2  in  the  queue  to  within 
K2  p2 


(l-hP2)  ■ 

Since  in  this  model  the  completion  times  of  the  type  2  customers  do 
not  form  a  sequence  of  independent  random  variables,  it  is  difficult  to 
say  much  about  their  expected  waiting  time  by  the  approaches  used  in  this 
paper. 
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Chapter  6 

MORE  GENERAL  QUEUES 


1 .  1 n  troduct ion 

V/e  now  relax 
wh i le  ma inta in  ing 
case  of  dependent 
expec  ted  wa i t  and 


some  of  the  independence  assumptions  of  previous  chapters 
stationar  ity.  Tandem  queues  are  taken  as  a  particular 
input.  In  this  case  an  expression  is  found  for  the  total 
is  seen  to  depend  on  some  unknown  covariance  terms. 


2 .  Stationary  Queues 

When  E[U]  <  0  and  the  sequence  {Un}  is  strictly  stationary  and  metrically 
transitive,  Loynes  (1962  (a))  has  shown  that  the  queue  has  a  stationary 
waiting  time  distribution.  Also  E[W]  and  E[ \~]  exist  and  equation  (1)  still 
holds  for  this  class  of  queues;  that  is, 


aoE[l]  =  -E[U3. 

Note  that  the  proof  of  theorem  1.1  required  stationar ity  but  not 
independence . 


3 .  Queues  with  Independent  Services 

If  the  service  times  form  an  independent  sequence  and  the  inter-arrival 
times  form  a  stationary  sequence,  from  equation  (2)  we  can  obtain  a  general¬ 
ization  of  theorem  1.2.  Using  the  same  method  of  proof,  but  realizing  that 

now  W  and  T  are  not  necessarily  independent  it  is  easy  to  show  that 
n  n 


E[W] 


EfU2!  Efl2! 
-2E[U]  "  2E[l] 


Cov  (W  ,T  ) 
x  n  n/ 

E[U] 


(52.) 
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Given  the  form  of  the  dependence  of  t he  inter-arrival  stream  it  may 
in  a  given  case  be  possible  to  calculate  the  covariance  term. 


4 .  Tandem  Queues 

Let  us  consider  now  m  single  channel  queues  in  tandem,  where  the  output 
from  one  becomes  the  input  into  the  next.  We  shall  assume  unlimited  queueing 
space  before  each  facility  so  that  no  blocking  occurs,  that  the  original 
input  forms  a  renewal  process,  and  that  the  service  times  in  each  facility 
are  independent  sequences  of  independent  random  variables.  That  is,  if 

5.  .  is  the  service  time  of  customer  l<  at  stage  i,  S  .  ~G.,  i=l,  2,  . . .  ,  m, 

K  y  I  K  j  I  I 

and  [S^  j}  are  mutually  independent. 

The  input  stream  to  a  1 1  facilities  after  the  first  will  not  in  general 

be  a  sequence  of  independent  random  variables  (the  case  of  Poisson  arrivals 

with  exponential  service  at  each  stage  is  the  exception).  Let 

2  .  .  . 
ct  =  variance  of  input  to  stage  i 

i 


a  =  variance  of  service  at  stage 
9i 


jo, .  =  service  rate  at  stage  i,  p.  = 


/  n  n 

'  =  n-th  idle  moment  at  stage  i 


i 

X_ 

P'i 


i 


w 


(i) 


T 

T 


n 

(i) 

n 

0) 

n 


=  waiting  time  of  n-th  customer  at  stage  i  (in  queue) 
=  inter-arrival  time  at  stage  ? 

=  inter-output  time  at  stage  i 


From  (  52)  we  immediately  obtain 


Note  that  there  is  cancellation  of  the  idle  time  moments  of  all  except 
the  last  facility  (this  happens  only  when  p  is  the  same  for  each  facility). 

The  completely  unknown  quantities  in  the  equation  are  the  covariance 
terms.  The  author  has  had  little  success  in  determining  their  order  of 
magnitude.  It  may  be  possible  to  find  bounds  on  these  quantities  in  some 
generality  since  each  term  is  the  covariance  between  the  wait  in  queue  of 
the  n-th  customer  and  the  time  until  the  (n+1 ) -th  customer  arrives.  it 
seems  to  the  author  that  in  general  this  correlation  would  be  small. 


summary 


The  aim  of  this  paper  is  to  find  simple  expressions  which  approximate 
some  of  the  measures  of  performance  in  the  Gi/G/1  queue.  The  large  body  of 
queueing  literature  shows  that  exact  expressions  for  many  of  these  measures 
are  extremely  complicated.  Often  they  are  implicit  in  nature,  making  them 
impractical  for  direct  application.  Notable  exceptions  are  the  papers  of 
Kingman  (1962  (a)  and  (b))  and  Newell  (1965).  In  one  paper  Kingman  deals 
with  the  GI/G/1  queue,  and  in  the  other  he  deals  with  asymptotic  properties 
(as  p  *-♦  1“)  of  stationary  queues.  Newell's  paper  deals  mainly  with  traffic 
light  problems. 

In  Chapter  1,  the  moments  of  the  waiting  time  and  idle  time  are  related 
by  equating  input  with  output  in  a  stationary  queue.  New  expressions  are 
found  for  the  mean  and  variance  of  the  wait,  and  an  expression  for  the 
mean  queue  length  follows  immediately.  Upper  bounds  for  all  Gl/G/1  queues 
are  found  easily  from  the  non-negativi ty  of  the  idle  time  variance  and  the 
bound  of  1  on  the  probability  a  customer  finds  the  system  empty  (Kingman 
(1962  (b))  finds  the  same  upper  bound  for  the  expected  wait).  These  bounds 
are  in  terms  of  the  means  and  variances  of  the  arrival  and  service  streams. 

A  non-trivial  lower  bound  on  the  expected  wait  is  found  which  requires 
knowledge  of  the  arrival  and  service  distributions. 

In  Chapter  2  we  recognize  that  an  idle  time  distribution  is  some 

complicated  tail  distribution  of  an  inter-arrival  time.  Restrictions  are 

placed  on  the  arrival  distribution  which  enable  us  to  obtain  "good"  bounds 

on  such  measures  as  the  mean  wait  and  mean  number  in  queue.  When  the  mean 

residual  life  of  an  inter-arrival  time  is  assumed  bounded  above  by  an 

r»  /\C  (u)  du  1 

ord  i na  ry  mean  i  nter-a  r r  i va  1  time,  that  is,  — ‘ ~  r-  a  1 1  T  =?  0  ,  we 
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find  bounds  which  give  the  mean  queue  length  to  within  When  the 

stronger  assumption  is  made  that  the  arrival  distribution  has  increasing 
failure  rate  (see  Chapter  2)  the  mean  queue  length  is  bounded  to  within 

ca+p  2 

( — 2~)  >  with  C1  -  i  (for  the  D/G/l  queue  this  reduces  to  “)  .  These  bounds 


are  in  terms  of  on  1  y  the  means  and  variances  of  the  arrival  and  service 
distributions.  Bounds  on  the  mean  idle  time  give  bounds  on  the  mean  length 
of  a  busy  period,  and  hence  on  the  mean  number  served  in  a  busy  period. 
Bounds  on  the  mean  actua  1  wait  give  bounds  on  the  mean  v i r tua 1  wa  j t .  Upper 
and  lower  bounds  are  also  found  when  the  mean  residual  life  of  an  inter- 


co  c  / 

arrival  times  is  bounded  below  (that  is,  \  - — g  a  1 1  T  ~  0)  , 

I  AC(T)  X 


and 


when  the  arrival  distribution  has  decreasing  failure  rate. 

Chapter  3  deals  with  these  generalizations  of  the  GI/G/1  queue.  Batch 
arrivals  are  treated  by  redefining  service  times,  and  when  A(t)  has  increasing 
failure  rate  or  has  mean  residual  life  bounded  above  by  the  mean  queue 
length  is  bounded  to  within  at  most  an  average  batch  size.  Queues  with 
service  in  batches  of  fixed  size  N  are  treated  by  redefining  inter-arrival 
times.  For  IFR  arrivals  the  mean  queue  length  is  bounded  to  within 


N 

2 


Queues  are  also  considered  where  the  service  breaks  down  in 


busy  periods  in  a  Poisson  manner. 

Queues  where  the  first  customer  in  each  busy  period  has  some  added 
delay  are  dealt  with  in  Chapter  4.  The  GI^/G/1  queue  (where  m  customers 
start  a  busy  period)  is  used  to  illustrate  the  results.  For  m  >  1  simple 
bounds  are  found  which  give  the  mean  queue  length  to  within  ~  when  A(t) 
is  IFR. 

Some  two-priority  queues  are  dealt  with  in  Chapter  5,  and  two  models 
are  considered.  In  the  first  one  the  arrivals  are  assumed  to  generate  two 


v-,v.  rir-' 
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renewal  processes,  and  high  priority  customers  are  Poisson.  In  the  second 
one  the  arrivals  are  assumed  to  generate  a  single  renewal  process,  and  the 
probability  that  an  arrival  is  a  high  priority  customer  is  p  ((1-p)  that 
he  is  a  low  priority  customer). 

The  last  chapter  deals  with  more  general  queues.  Stationarity  is 
retained  but  the  independence  assumptions  are  weakened.  Some  preliminary 
results  on  tandem  queues  are  given,  but  the  expression  for  the  expected 
wait  has  some  covariance  terms  of  unknown  order  of  magnitude. 

There  is  much  to  be  done  in  the  are  of  approximations  in  queueing. 

The  importance  of  the  idle  time  distribution  is  clearly  demonstrated  in 
this  paper.  More  work  needs  to  be  done  in  relating  its  properties  to  those 
of  the  arrival  and  service  streams. 
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